Walking into an Absolute Sum 
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1 Introduction 

Recently, it was asked by Paul Bruckman 1 to show that the sum 



%\ S r (n) = J2(^)\n-k\ r (1) 



2n / n 

In 



k=0 



evaluates to n 2 i^^) ^ or r = 3. In the published solution [TJj|, it was also noted that S\[n) = n ( 2 ^J j 

and, as a consequence, it was conjectured that S 2r+ i(n) equals the product of y^j and a monic 
polynomial of degree r + 1. 

We show this conjecture to be true, albeit with the modification of discarding the adjectival 
' modifier "monic." In fact, we show that 

g '. S 2 r+i(n) = P r (n)n( 2 ™) and S 2r (n) = Q r (n) 2 2n ~ r , 

where P r {n) and Q r (n) are both polynomials of degree r with integer coefficients. We then in- 
vestigate the relationship of these polynomials to the Dumont-Foata polynomials [5]. These are 
generalizations of the Gandhi polynomials, which find their origin in a representation of the Genoc- 
chi numbers, first conjectured by Gandhi [9]. Finally, we show that the sums S r (n) are essentially 
the moments of a random variate, measuring the absolute distance to the origin in a symmetric 
Bernoulli random walk, after In time steps. 



2 Derivation 

We note that the sum can be rewritten as 



rt 



^)= 2 E( n 2 :> r -( 2 ;)^ 



n 



k=0 

with c) r o the Kronecker delta. Now consider, for r > 1, 

n'S r (n)-S r+2 (n) = 2j2( ^ ,) k r U - k 2 ) = 4n(2n - 1) ]rV \ V' 

leading directly to the recursion 

S r+2 (n) = n 2 S r (n) - 2n(2n - l)5 r (n - 1). (2) 
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For r = the derivation is slightly more elaborate, because we need to keep track of the additional 
term, but leads to the same recursion, so that © is valid for all nonnegative integers r. To start 
the recursion, we find the value So(n) = 2 2n by an application of the binomial theorem to The 
value of Si(n) is easily obtained by breaking up the summand k to create two sums: 



Si(n) = J2 

k=0 



2n 
n — k 



(n + k) — (n — k) 



2n£ 

fc=0 



2n 
n - 
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n-X 
k=0 



In - 1 
n — k — 1 



and one sees that, after changing the range of summation of the second sum to start at k = 1, 
all terms cancel out, with the exception of the summand 2n ( 2n n ~ 1 ^- Rearranging terms gives the 

desired Si(n) = n ( 2 ™)- 

It is now clear that the structure of the sum depends upon the parity of r. Starting with the 
odd values, we simplify the recursion © by the substitution 5 , 2 r +i(n) = P r (n)n ( 2 ^) to give 



P 



■+i(») = 



= n 2 



P r {n) -P r (n- 1)1 + nP r {n - 1), 



(3) 



with initial condition Po(n) = 1. An inductive argument now shows that P r (n) is a polynomial of 
degree r with integer coefficients, and proves the modified conjecture. It is not difficult to show 
that r! is the leading coefficient of P r (n), and, hence, that these polynomials are not monic. In fact, 
the only cases for which the leading coefficient is 1 are r = and r = 1. The first few polynomials 
are now easily determined as: 

P (n) = 1, 

P\{n) = n, 

P 2 (n) = (2n - l)n, 

P 3 (n) = (6n 2 -8n + 3)n, 

P 4 (n) = (24n 3 - 60n 2 + 54n - 17)n, 

P 5 (n) = (120n 4 - 480n 3 + 762n 2 - 556n + 155)n. 

For the even sums we substitute S2r(n) = Q r (n) 2 2n ~ r to give the recursion 

Q r+ x{n) = 2n 2 Q r {n) - Q r (n- 1) +nQ r (n-l), (4) 

with initial condition Qo(n) = 1. This shows that Q r (n) is a polynomial of degree r with integer 
coefficients. It is not difficult to establish that the leading coefficient is given by (2r — 1) • (2r — 
3) • • • 3 • 1 = (2r)!/(2 r r!), and, hence, that these polynomials are also not monic. Applying the 
recursion gives the first few polynomials as: 



Qo(n 
Qi(n 
Q2(n 
Q?,{n 

Q 5 (n 



1, 

n, 

(3n — l)n, 
(15n 2 - 15n + 4)n, 
(105n 3 - 210n 2 + 147n 
(945n 4 



3150n 3 + 4095n 2 



34)n, 

2370n + 496)n. 



It is worth noting that, by evaluating S r (n) for particular values of n, one can derive various 
properties of [the coefficients of] the polynomials P r (n) and Q r (n). For instance, it is not difficult 
to show that the coefficients of P r {n) sum to unity, and those of Q T {n) to 2 r ~ 1 (for r > 1), by 
evaluating the sums for n = 1. Indeed, one can derive the closed- form solutions for S2 r (n) and 
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S 2r +i(n), by solving a system of linear equations in r unknowns, representing the coefficients of 
the corresponding polynomial. 

In the constant of the polynomials P r {n)/n one recognizes the Genocchi numbers (UE3) named 
after the Italian mathematician Angelo Genocchi (1817-1889): 

G 2 = -l, G 4 = l, G 6 = -3, G 8 = 17, Gi = -155, G 12 = 2073, .... 

These integers are defined through the exponential generating function 

2t t 2r 

and are related to the Bernoulli numbers by G 2r = 2(1 — 2 2r )B 2r . The Genocchi numbers are 
listed as sequence A001469 in the on-line version of the encyclopedia of integer sequences |15j . 
where additional references may be found. The constant of the polynomials Q r (n)/n matches the 
first terms of the sequence A002105 in which is generated by 2 r ~ 1 G 2r /r, and related to the 
tangent numbers. The connection to the Genocchi numbers will be further explored in the next 
section, where the polynomials P r (n) and Q r (n) are found to be related to special cases of the 
Dumont-Foata polynomials. 

Another matter of interest is the leading coefficient of the polynomials, characterizing the 
behavior of the sums S r (n) for large values of n. For the even indexed sums, this is easily established 
as 

^(n)~|^2 2 ^, (5) 
and for the odd indexed sums we can use Stirling's formula to give ( 2 ^) ~ 2 2n / yjim, so that 

S 2r+1 (n)~^=2 2 "n r '+i (6) 
V 71 " 

In these expressions, one can recognize the moments of a central chi-distribution, see for instance |12l 
pp. 420-421]. That this is no coincidence will be shown in section 0] where we establish the con- 
nection between the sums S r {n) and the distance to the origin in a symmetric Bernoulli random 
walk. 



3 Dumont-Foata polynomials 

In this section we show that the polynomials P r (n) and Q r (n) are related to special cases of the 
Dumont-Foata polynomials [H]. These are defined recursively by means of 

F r+1 (x,y,z) = (x + z){y + z)F r (x,y, z + 1) - z 2 F r (x,y,z), 

with initial condition F\ (x, y,z) = 1. Explicit expressions for these polynomials and their generating 
functions have been derived by Carlitz [2], but are too lengthy to display here. 

The Dumont-Foata polynomials can be regarded as generalizations of the Gandhi polynomials, 
see for instance jSJ El > which are defined by the recursion 

P r+ l{z) = {Z + l) 2 P r (z + 1) - Z 2 P r {z), 

with initial condition P\{z) = 1. The coefficients of the first few of these polynomials are displayed 
in Tabled and can also be found in [151 Seq. A036970]. The Gandhi polynomials arose from a 
conjecture made by Gandhi 9 , concerning a representation of the Genocchi numbers. Gandhi's 
conjecture that P r (0) = (—l) r G 2r was proved by Carlitz 0, and also by Riordan and Stein [Ti] . 
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1 

2 1 
6 8 3 

24 60 54 17 

120 480 762 556 155 

720 4200 10248 12840 8146 2073 



Table 1: Coefficients of the Gandhi polynomials, arranged in triangular form. 

1 

3 1 
15 15 4 

105 210 147 34 

945 3150 4095 2370 496 

10395 51975 107415 111705 56958 11056 

Table 2: Coefficients of the polynomials Q r (z), arranged in triangular form. 

Another polynomial that can be derived as a special case of the Dumont-Foata polynomials is given 
by Q r (z) = 2 r '~ 1 F r (i, 1, z), and is generated by the recursion 

Q r+ l(z) = (2Z + 1)0 + 1)4.(2 + 1) - 2z 2 Qr{z), 

with initial condition Q\{z) = 1. The coefficients of the first few of these polynomials are 
displayed in Table I2J For r > 1, one can easily verify by substitution in (J3J) and ® that 
P r (n) = (—l) r ~ 1 nP r (—n) and Q r (n) = (— l) r ~ l nQ r {— n). This gives the connection to the Dumont- 
Foata polynomials (for positive r) as: 

P r (n) = (-l) p_1 nF r (l,l,-n) 

and 

Q r (n) = (-2) p - 1 n J F r (^,l J -n). 

The occurrence of the Genocchi numbers in the constant of the polynomials P r (n)/n is now seen 
to be a direct consequence of Gandhi's conjecture that F r (l,l,0) = (—VfG<2 r . The occurrence of 
the Genocchi numbers in the constant of the polynomials Q r {n)/n is conjectured by the present 
author in the form F r (^, 1,0) = (— l) r G2 r /r. 

4 Symmetric Bernoulli random walks 

In a symmetric Bernoulli random walk, one considers the movements of a particle starting at time 
t = at the origin. Its movements are determined by a chance mechanism, where a fair coin is 
flipped and the particle is moved one unit to the right if it is heads up, and one unit to the left 
if it is tails up. A more exhaustive description and in-depth study of random walks can be found 
in Feller |5| or Revesz |13| . A more playful introduction to the topic is given in the monograph by 
Dynkin and Uspenskii A topic of interest is the position of the particle after 2n coin tosses: 
Ym = X\ + X2 + • • • + A^2n, where Xi is +1 or — 1 depending upon whether or not the coin showed 
heads in the ith coin toss. Note that the Xi are independent and identically distributed variates 
with mean and variance 1. The probability distribution of the position of the particle after 2n 
moves can be derived from a simple combinatorial argument, see for instance [HI p. 75] or |13| p. 13], 
and is given by 

2n N ~-2n 



Prob{Y 2n = 2k) = [ n _ k )2 
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where k = —n, —n + 1, . . . , n, and n a positive integer. The matter of interest in the context of this 
note is the distance to the origin, |^2n| at time t = 2n. Its moments are given by 

and one sees that E|Y2 n | r = 2 r ~ 2n S r (n), thus establishing the connection to the absolute sums 
from the introduction. The limit behavior of these sums now becomes clear. By the central limit 
theorem, see for instance [111 p. 18], one has that Yi n , for sufficiently large n, follows a normal 
distribution with mean and variance 2n. This implies that asymptotically, |^2 n | has a half- normal 
or central chi-distribution, so that 

I m\ T(l/2) 

see for instance [121 pp. 420-421]. This gives the asymptotic behavior of the sums as 

S r (n) = 2 2n ~ r E \Y 2n \ r ~ r[(r / + / 1) / 2] 2 2n n r / 2 , 

rv J I zn\ r(l/2) 

and upon expanding the gamma functions one recovers the limit results (jSJ) and © . 

5 Discussion 

One could possibly use the relation of the Gandhi polynomials to the sums SW+i^n,) to gain new 
insight into the former. In particular, one now has an expression to derive the function values of 
the Gandhi polynomials for negative, integral arguments: 

aho = ( 2 ;) - 1 ± (;_-,) *~>. 

For example, one easily obtains P r (-1) = (-l) r_1 and P r (-2) = (-l) r ^ 1 (2 2r ^ 1 + l)/3. 

Likewise, one can use the relation of the moments of the absolute distance to the origin in a 
symmetric Bernoulli random walk and the sums S r (n) to express these moments in terms of the 
polynomials P r (n) and Q r (n): 



E\Y 2n \ 2r = 2 r Q r (n) and E \Y 2n \ 2r+1 = ( ln ) 2 2 ^ +1 n P r (n). 



(In 
V n 

This equivalence can be used to establish the rate of convergence to the moments of the half-normal 
distribution. 

Finally, it should be noted that one can also determine expressions for S2 r (n) by means of the 
generating function 



2n / n 

2n 



fM = E 7 ein ~ h)v = e ^ t 1 + e ~T n = 2 n [i + cosier , 

so that Szrin) = fn r \o). However, this approach covers only the even indexed case, and does not 
give the same insight into the problem as the one that we have followed here. 
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